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Mycobacterium tuberculosis (Mtb), the causative agent of 
human tuberculosis, remains one of the most prevalent 
human pathogens and a major cause of mortality world- 
wide. Metabolic network is a central mediator and defin- 
ing feature of the pathogenicity of Mtb. Increasing evi- 
dence suggests that lysine succinylation dynamically 
regulates enzymes in carbon metabolism in both bacteria 
and human cells; however, its extent and function in Mtb 
remain unexplored. Here, we performed a global succi- 
nylome analysis of the virulent Mtb strain H37Rv by using 
high accuracy nano-LC-MS/MS in combination with the 
enrichment of succinylated peptides from digested cell 
lysates and subsequent peptide identification. In total, 
1545 lysine succinylation sites on 626 proteins were iden- 
tified in this pathogen. The identified succinylated pro- 
teins are involved in various biological processes and a 
large proportion of the succinylation sites are present on 
proteins in the central metabolism pathway. Site-specific 
mutations showed that succinylation is a negative regu- 
latory modification on the enzymatic activity of acetyl- 
CoA synthetase. Molecular dynamics simulations demon- 
strated that succinylation affects the conformational 
stability of acetyl-CoA synthetase, which is critical for its 
enzymatic activity. Further functional studies showed that 
CobB, a sirtuin-like deacetylase in Mtb, functions as a 
desuccinylase of acetyl-CoA synthetase in in vitro assays. 
Together, our findings reveal widespread roles for lysine 
succinylation in regulating metabolism and diverse pro- 
cesses in Mtb. Our data provide a rich resource for 
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Post-translational modifications (PTMs)' are complex and 
fundamental mechanisms modulating diverse protein proper- 
ties and functions, and have been associated with almost all 
known cellular pathways and disease processes (1, 2). Among 
the hundreds of different PTMs, acylations at lysine residues, 
such as acetylation (3-6), malonylation (7, 8), crotonylation (9, 
10), propionylation (11-13), butyrylation (11, 13), and succi- 
nylation (7, 14-16) are crucial for functional regulations of 
many prokaryotic and eukaryotic proteins. Because these 
lysine PTMs depend on the acyl-CoA metabolic intermedi- 
ates, suchas acetyl-CoA (Ac-CoA), succinyl-CoA, and malonyl- 
CoA, lysine acylation could provide a mechanism to respond 
to changes in the energy status of the cell and regulate energy 
metabolism and the key metabolic pathways in diverse or- 
ganisms (17, 18). 

Among these lysine PTMs, lysine succinylation is a highly 
dynamic and regulated PTM defined as transfer of a succinyl 
group (-CO-CH;-CH;-CO-) to a lysine residue of a protein 
molecule (8). It was recently identified and comprehensively 
validated in both bacterial and mammalian cells (8, 14, 16). It 
was also identified in core histones, suggesting that lysine 
succinylation may regulate the functions of histones and af- 
fect chromatin structure and gene expression (7). Accumulat- 
ing evidence suggests that lysine succinylation is a 
widespread and important PTM in both eukaryotes and pro- 
karyotes and regulates diverse cellular processes (16). The 
system-wide studies involving lysine-succinylated peptide 


! The abbreviations used are: TB, Tuberculosis; Mtb, Mycobacte- 
rium tuberculosis; PTM, post-translational modification; Acs, acetyl- 
CoA synthetase; FDR, false discovery rate; GO, Gene Ontology; 
KEGG, Kyoto Encyclopedia of Genes and Genomes; SPC, simple 
point charge; LINCS, linear constraint solver; PME, particle mesh 
Ewald; MD, molecular dynamics; Ac-CoA, acetyl-CoA; CoA, Coen- 
zyme A; NAD', nicotinamide adenine dinucleotide; NAM, Nicotin- 
amide; Sirt5, sirtuin; RMSD, root mean square deviation; Rgs, gyra- 
tions; RMSF, residue-specific root mean square fluctuation. 
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immunoprecipitation and liquid chromatography-mass spec- 
trometry (LC-MS/MS) have been employed to analyze the 
bacteria (E. coli) (14, 16), yeast (S. cerevisiae), human (HeLa) 
cells, and mouse embryonic fibroblasts and liver cells (16, 19). 
These succinylome studies have generated large data sets of 
lysine-succinylated proteins in both eukaryotes and pro- 
karyotes and demonstrated the diverse cellular functions of 
this PTM. Notably, lysine succinylation is widespread among 
diverse mitochondrial metabolic enzymes that are involved in 
fatty acid metabolism, amino acid degradation, and the tri- 
carboxylic acid cycle (19, 20). Thus, lysine succinylation is 
reported as a functional PTM with the potential to impact 
mitochondrial metabolism and coordinate different metabolic 
pathways in human cells and bacteria (14, 19-22). 

Mycobacterium tuberculosis (Mtb), the causative agent of 
tuberculosis (TB), is a major cause of mortality worldwide and 
claims more human lives annually than any other bacterial 
pathogen (23). About one third of the world’s population is 
infected with Mtb, which leads to nearly 1.3 million deaths and 
8.6 million new cases of TB in 2012 worldwide (24). Mtb 
remains a major threat to global health, especially in the 
developing countries. Emergence of multidrug resistant 
(MDR) and extensively drug-resistant (XDR) Mtb, and also the 
emergence of co-infection between TB and HIV have further 
worsened the situation (25-27). Among bacterial pathogens, 
Mtb has a distinctive life cycle spanning different environ- 
ments and developmental stages (28). Especially, Mtb can 
exist in dormant or active states in the host, leading to asymp- 
tomatic latent TB infection or active TB disease (29). To 
achieve these different physiologic states, Mtb developed a 
mechanism to sense diverse signals from the host and to 
coordinately regulate multiple cellular processes and path- 
ways (30, 31). Mtb has evolved its metabolic network to both 
maintain and propagate its survival as a species within hu- 
mans (32-35). It is well accepted that metabolic network is a 
central mediator and defining feature of the pathogenicity of 
Mtb (23, 36-38). Knowledge of the regulation of metabolic 
pathways used by Mtb during infection is therefore important 
for understanding its pathogenicity, and can also guide the 
development of novel drug therapies (39). On the other hand, 
increasing evidence suggests that lysine succinylation dy- 
namically regulates enzymes in carbon metabolism in both 
bacteria and human cells (14, 19-22). It is tempting to spec- 
ulate that lysine succinylation may play an important regula- 
tory role in metabolic processes in Mib. However, to the best 
of our knowledge, no succinylated protein in Mtb has been 
identified, presenting a major obstacle to understand the reg- 
ulatory roles of lysine succinylation in this life-threatening 
pathogen. 

In order to fill this gap in our knowledge, we have initiated a 
systematic study of the identities and functional roles of 
the succinylated protein in Mtb. Because Mtb H37Rv is the 
first sequenced Mtb strain (40) and has been extensively used 
for studies in dissecting the roles of individual genes in patho- 


genesis (41), it was selected as a test case. We analyzed the 
succinylome of Mtb H37Rv by using high accuracy nano-LC- 
MS/MS in combination with the enrichment of succinylated 
peptides from digested cell lysates and subsequent peptide 
identification. In total, 1545 lysine succinylation sites on 626 
proteins were identified in this pathogen. The identified suc- 
cinylated proteins are involved in various biological processes 
and render particular enrichment to metabolic process. A 
large proportion of the succinylation sites are present on 
proteins in the central metabolism pathway. We further dis- 
sected the regulatory role of succinylation on acetyl-CoA 
synthetase (Acs) via site-specific mutagenesis analysis and 
molecular dynamics (MD) simulations showed that reversible 
lysine succinylation could inhibit the activity of Acs. Further 
functional studies showed that CobB, a sirtuin-like deacety- 
lase in Mtb, functions as a deacetylase and as a desuccinyl- 
ase of Acs in in vitro assays. Together, our findings provide 
significant insights into the range of functions regulated by 
lysine succinylation in Mtb. 


EXPERIMENTAL PROCEDURES 


Bacterial Strains and Growth Conditions—The M. tuberculosis 
strain H37Rv (ATCC 27294) was obtained from the National Center for 
Medical Culture Collections (CMCC, Beijing, China). The seed culture 
was grown to log phase (ODgo9 ~ 1.2) in Middlebrook 7H9 liquid 
medium (Difco) supplemented with ADC (Albumin-dextrose-catalase) 
at 37 °C. After 7 days of growth of log phase, the cultures were pellet, 
washed twice with PBS (0.1 M Na;HPO,, 0.15 M NaCl, pH 7.5), and 
resuspended in fresh liquid medium followed by incubation for 2-3 
weeks to a final ODg99 ~ 0.6-0.8. For carbon source utilization 
studies, Mtb seed culture was resuspended in fresh complex growth 
media containing additional carbon sources, such as 0.2% pyruvate, 
0.2% glucose, 0.2% glycerol, 0.2% succinate, or 0.2% acetate, 
respectively. 

Preparation of Protein Lysate and In-solution Trypsin Digestion—To 
identify as many succinylation events as possible, we have performed 
a series of preliminary experiments using different amount of Mtb 
proteins to find the best protocol. Following is the optimized experi- 
mental procedures based on these preliminary experiments. In brief, 
bacterial cells (250 ml) were harvested by centrifugation at 2500 x g 
for 10 min at 4 °C. The pellet was washed twice with ice-cold PBS, 
resuspended in ice-cold PBS containing 1 protease inhibitor mix- 
ture (Roche Diagnostics Ltd, Mannheim, Germany) and 1 mw PMSF 
(Beyotime Institute of Biotechnology, Jiangsu, China). The mixture 
was placed in Lysing Matrix B bead beater vials (pre-filled with 0.1 
mm silica; MP Biomedicals, Santa Ana, CA) and homogenized using 
the FastPrep-24 sample preparation system (MP Biomedicals, CA). 
The cells were lysed by bead-beating for five cycles of 35 s at a 
setting of 6.5 m/s. Between each cycle of bead beating, the samples 
were incubated on ice for 3 min. Cell lysates were cleared of silica and 
cellular debris via centrifugation at 12,000 X g for 15 min at 4 °C. 
Supernatant was then filtered through the 0.2 um syringe filters twice 
to sterilize lysate. The clarified lysate was precipitated by using 1096 
trichloroacetic acid and 196 sodium deoxycholate, and then washed 
twice with ice cold acetone. The precipitated proteins (2 mg) were 
redissolved in 50 mm ammonium bicarbonate, then in-solution di- 
gested by trypsin according to previously described (14, 42). The 
protein concentration was determined with the Bradford assay (Bio- 
Rad, Hercules, CA). 
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Immunoaffinity Enrichment of Lysine Succinylated Peptides — Suc- 
cinylated peptides were enriched using agarose-conjugated anti-suc- 
cinyl lysine antibody (PTM Biolabs Inc., Chicago, IL) as previously 
described (14). Briefly, tryptic peptides were redissolved in NETN 
buffer (50 mm Tris-HCI, pH 8.0, 100 mm NaCl, 1 mm EDTA, 0.5% 
Nonidet P-40), and then incubated with anti-succinyllysine antibody 
conjugated protein A agarose beads at 4 °C for 6h with gentle rota- 
tion. The supernatant was removed and the beads were washed three 
times with NETN buffer, twice with ETN buffer (50 mm Tris-Cl, pH 8.0, 
100 mm NaCl, 1 mM EDTA), followed by three times wash with water. 
The bound peptides were eluted by washing three times with 196 
trifluoroacetic acid. The eluates were combined and dried in a Speed- 
Vac. The resulting succinylated peptides were loaded onto self- 
packed C44 STAGE tips according to the manufacturer's instructions 
to desalt the sample, prior to nano-HPLC-MS/MS analysis. 

Nano-HPLC-MS/MS Analysis— The enriched peptides were dis- 
solved in the HPLC buffer A (0.196 (v/v) formic acid in water), and 
analyzed by online nanoflow LC-MS/MS using an easy nLC-1000 
system (Thermo Scientific) connected to a Q-Exactive (Thermo Sci- 
entific) mass spectrometer. Briefly, samples were injected onto the 
analytical C,,-nanocapillary LC column (5 um particle size, 100 A 
pore diameter) and eluted at a flow rate of 300 nl/min with a 45 min 
gradient from 5% solvent B (9096 ACN/0.1% formic acid, v/v) to 80% 
Solvent B. The peptides were then directly ionized and sprayed into a 
Q-Exactive mass spectrometer by a nanospray ion source. Mass 
Spectrometer was operated in a data-dependent mode with an auto- 
matic switch between MS and MS/MS acquisition. Full MS spectra 
from m/z 350 to 1600 were acquired with a resolution of 70,000 at 
m/z — 200 in profile mode. Following every survey scan, up to 15 
most intense precursor ions were picked for MS/MS fragmentation by 
higher energy C-trap dissociation (HCD) with normalized collision 
energy of 3596. Lock mass at m/z 445.12003 was enabled for full MS 
scan. The dynamic exclusion duration was set to be 5 s with a repeat 
count of one and + 10 ppm exclusion window. 

Protein Sequence Database Search — All acquired raw data were 
processed with the MaxQuant software (version 1.3.0.5) (43). The 
peak lists were searched against the Mtb H37Rv protein database 
from NCBI ftp site (http://www.ncbi.nlm.nih.gov/) containing 4018 
protein sequences. Two missed cleavages were allowed for trypsin. 
The precursor and fragment ion mass tolerances were 10 ppm and 
0.02 Da, respectively. Carbamidomethylation (Cys) was set as a fixed 
modification, whereas oxidation (Met), deamidation (Asn/Gln), succi- 
nylation (Lys), and acetylayion (Protein N-terminal) were set as vari- 
able modifications. Minimum peptide length was set at six. The esti- 
mated false discovery rate (FDR) thresholds for modification site, 
peptide and protein were specified at maximum 196. All MS/MS 
spectra for the identified peptides with succinylation modifications 
were manually inspected using criteria as previously reported (44, 45). 
Furthermore, to improve the reliability of result, all succinylation sites 
assigned to the peptide C-terminal were removed, prior to bioinfor- 
matics analysis. 

Bioinformatics Analysis — The identified succinylated proteins were 
grouped into biological process and molecular function class based 
on the Gene Ontology (GO) terms by Blast2GO software (46). The 
subcellular localization of succinylated proteins was analyzed with 
PSORTb program (47). To gain information on the over-representa- 
tions, GO enrichment analysis was performed using Cytoscape plugin 
BiNGO (48, 49) with the parameters as previously described (50). The 
reference GO ontology in Cytoscape ontology format was created 
using whole Mtb H37Rv proteome that extracted from MTBbase 
(http://www.ark.in-berlin.de/Site/MTBbase.html). Pathway and Pfam 
functional domain enrichment analyses were performed using the 
DAVID bioinformatics resource (51, 52). Functional interaction net- 
work analysis was performed using interaction data as previously 


reported (53) and the interaction network was visualized by Cyto- 
scape (version 3.0.2) (48). Amino acid sequence motifs were analyzed 
using motif-X (64) and an in-house script as previously reported (42). 
The position-specific heatmap was generated by plotting the log,, of 
the ratio of frequencies. The secondary structure types were deter- 
mined for positions with succinylated lysine using NetSurfP tool (55). 
In our data, the corresponding p value < 0.05 was considered sta- 
tistically significant. 

Production of Specific Antibodies against Mtb Proteins—The gen- 
erations of polyclonal antibodies against the specific Mtb proteins 
were carried out by ABclonal Inc. (Wuhan, Hubei, China). Briefly, to 
produce antibodies against the Mtb proteins, i.e. ATP synthase alpha 
chain (AtpA), ATP synthase beta chain (AtpF), 3-oxoacyl-[acyl-carrier 
protein] synthase 1 (KasA), acyl-[acyl-carrier protein] desaturase 
(DesA1), acyl-CoA ligase (FadD31), F420-dependent glucose-6- 
phosphate dehydrogenase (Fgd1), malate synthase G (GlcB), heat 
shock protein (HspX), two component system transcriptional regula- 
tor (MtrA) and polyketide synthase (Pks2), the full-length of the genes 
were cloned into the expression vector pGEX-4T (Pharmacia, Pisca- 
taway, NJ), and then each expression vector was transformed into 
E. coli strain BL21 (DES) for overexpression. For expression, 1 mM 
isopropyl-8-b-1-thiogalactopyranoside (IPTG) was added into the 
logarithmically growing bacterial culture and cells were harvested 
after 4 h of induction at 30 °C. Following purification of these anti- 
gens, immunization and sampling of the anti-sera from rabbit were 
performed by ABclonal Inc. The specificity of the generated antibod- 
ies was determined by the manufacturer using ELISA and Western 
blotting. 

Immunoprecipitation and Western blotting—To validate succinyla- 
tion of identified proteins, we performed immunoprecipitation (IP) 
experiments to analysis the succinylated proteins in vivo using the 
generated specific protein antibodies. For IP, the protein specific 
antibody was conjugated to the protein G Dynabeads (Invitrogen AS, 
Oslo, Norway) in PBS by gentle rocking for 4 h at 4 °C. The conju- 
gated beads were washed three times with PBS and incubated with 
the whole cell lysates overnight at 4 °C. The beads were then washed 
three times with PBS to remove the unbound proteins. Bound pro- 
teins were boiled in SDS loading buffer for 5 min, then subjected to 
1596 SDS-PAGE and transferred to a polyvinylidene difluoride (PVDF) 
membrane. For Western blotting, the membrane was blocked over- 
night at 4 °C or 37 °C for 2 h in TBS buffer (25 mM Tris-HCl, pH 8.0, 
150 mM NaCl) containing 5% bovine serum albumin (BSA) and incu- 
bated with either the specific Mtb protein antibody or the anti-succinyl 
lysine antibody (PTM Biolabs Inc., Chicago, IL) (1:2000, in TBS/596 
BSA) overnight at 4 °C. After washing three times with TBST buffer 
(25 mm Tris-HCl, pH 8.0, 150 mm NaCl, 0.196 Tween20), the mem- 
brane was incubated with horseradish peroxidase- conjugated goat 
anti-rabbit antibody (1:5000 dilutions) for 1 h at 37 °C. The mem- 
brane was then washed with TBST buffer and visualized with en- 
hanced chemiluminescence (ECL) immunoblotting detection re- 
agents (Advansta Inc., Menlo Park, CA). The density of each band 
was determined with a fluorescence scanner (ImageQuant TL, GE 
Healthcare). Densitometry analysis was performed via the ImageJ 
suite (http://rsbweb.nih.gov/ij/). 

Site-directed Mutagenesis and Purification of Acs— The deacety- 
lase (CobB) and Acs expression vectors (pET32a-cobB and pET20b- 
acs, kindly provided by Professor Jiao-Yu Deng, Wuhan Institute of 
Virology, CAS) were transformed into E. coli BL21 (DE3) for protein 
expression. The site-directed mutation of acs gene was introduced 
into the selected sites by PCR reaction. The PCR product was trans- 
formed into E. coli BL21 (DES), and subsequently the mutagenesis 
was completely sequenced to confirm the presence of the site-di- 
rected mutation. Mutagenic primers were given below with the 
mutated base triplets underlined: K193R-sense: 5'-CGCGGCAGAC- 
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CATCGCCCCTCA AGGCGGCC-3', K193R-antisense: 5'-CGATGG- 
TCTGCCGCGCCGAAACTGCCCGTC-3', K366R-sense: 5'-TTT- 
ATGAGATGGGGCCGTGAGATCCCCGAC-3’, K366R-antisense: 5'- 
GCCCCATCTCATAAACATCCGGATGAGGGT-3’. The wild type 
E. coli BL21 (DES)/pET32a-cobB, E. coli BL21 (DE3)/pET20b-acs and 
the site-directed mutants of acs gene were grown in 5 ml Luria- 
Bertani (LB) medium containing ampicillin (50 ug/ml) overnight at 
37 °C, and then the cultures were transformed into 500 ml fresh LB 
medium with ampicillin (50 ug/ml) at 37 °C in shaking flasks to optical 
density at 600 nm of 0.4~0.6. Cells were induced with 0.2 mM of IPTG 
at 18 °C for 12h. The cultured cells were harvested by centrifugation 
and washed twice with ice-cold PBS, resuspended in binding buffer 
(20 mM Tris-HCl, pH 7.9, 500 mm NaCl and 5 mm imidazole). The 
mixture was disrupted by high pressure cell disrupter with an output 
of 1500 W (1500 bar, JN-02C, JNBIO, Guangzhou, China). Cellular 
debris was removed by centrifugation at 8000 x g for 30 min at 4 °C, 
and the resulting supernatants were loaded onto an affinity Ni?* 
column (Qiagen Inc., Chatsworth, CA) pre-equilibrated with the bind- 
ing buffer. After the column was washed with binding buffer, followed 
by washing buffer (20 mM Tris-HCl, pH 7.9, 500 mm NaCl and 60 mM 
imidazole or 80 mM imidazole), the protein was eluted with elution 
buffer (20 mM Tris-HCl, pH 7.9, 500 mm NaCl and 250 mm imidazole). 
The elution was collected and concentrated using a 30,000 MWCO 
concentrator (Millipore) in storage buffer (20 mm Tris-HCl, pH 8.0, 50 
mm NaCl). The protein concentration was determined with the Brad- 
ford assay (Bio-Rad). 

Acs Activity Assay—The activity of purified Acs and site-directed 
mutants of K193R, K366R, K193R/K366R proteins were measured as 
described (56, 57). Briefly, each reaction contained 50 mM Tris-HCl 
(pH 8.0), 600 mm hydroxylamine-HCI (pre-neutralized with KOH), 10 
mM MgCl, 2 mM pL-Dithiothreitol (DTT), 1 mm Coenzyme A (CoA), 30 
mM ATP, and 5 mM potassium acetate. Reactions were pre-incubated 
at 37 °C for 5 min before the addition of purified proteins (1 uM). 
Samples without protein served as a blank, and formation of acetyl- 
hydroxamate by acetyl-phosphate (AcP) served as a standard. After 
10 min incubation, the reactions were quenched by adding an equal 
volume of stop solution [2.5% (w/v) FeCl, in 2 M HCl, 10% (w/v) 
trichloroacetic acid]. The reaction products were centrifuged to re- 
move turbidity, and then measured at 540 nm. 

In Vitro Desuccinylation Activity Assay — The purified recombinant 
CobB and Acs were incubated in 50 yl of reaction buffer (50 mM 
Tris-HCl, pH 8.5, 187 mm NaCl, 2.7 mm KCI, 1 mm MgCl.) with or 
without 500 uM NAD+ and in the presence or absence of 2 mM 
Nicotinamide (NAM). The mixture was incubated at 25 °C for 10 h. 
The reaction products were analyzed by SDS-PAGE and Western 
blotting. 

Molecular Dynamics Simulations— The initial model of the Acs 
module was prepared from the available crystal structure (PDB code 
2P20, 50.5696 sequence identity) via structural modeling using CPH- 
models 3.2 Server (http://www.cbs.dtu.dk/services/CPHmodels). All 
simulations were based on the initial model and a total of seven 
systems were prepared for simulation using GROMACS 4.5.5 (58) in 
conjunction with the OPLS-AA/L all-atom force field (59). The protein 
was then solvated in simple point charge (SPC) water molecules in a 
cubic box (60), with the box edges —1.0 nm from any atom of the 
protein, and eighteen additional Na* ions were added to neutralize 
the charge of each system. Each system was then energy minimized 
using a steepest descents integrator either until the maximum force 
was less than 1000 kJ/mol/nm on any atom or until additional steps 
resulted in a potential energy change of less than 1 kJ/mol. Next, the 
simulations were performed under a constant temperature of 300 K 
and the V-rescale algorithm (61) was used with a temperature cou- 
pling time constant of 0.1 ps. All bond lengths were constrained using 
the linear constraint solver (LINCS) algorithm (62). Van der Waals 


interactions used a simple cut off at 1.4 nm and long-range electro- 
static interactions were handled using the particle mesh Ewald (PME) 
method (63, 64) with a fourth-order spline interpolation and a 0.1 nm 
Fourier grid spacing. Once each system was sufficiently equilibrated 
around the target temperature, isotropic was used for pressure cou- 
pling and the constant pressure was set to 1.0 bar in all directions and 
a pressure coupling time constant of 1.0 ps. Finally, each system was 
subjected to 10 ns of molecular dynamics (MD) simulation and the 
time step used in the simulations was 2 fs. All analyses were per- 
formed using the GROMACS suite of tools and a secondary structure 
recognition algorithm (DSSP) (65), which was implemented in 
GROMACS. The PyMOL Molecular Graphics System (version 1.7.2, 
http://www.pymol.org) was employed to present the structural results 
of this study. 


RESULTS 


Detection of Protein Lysine Succinylation in Mtb — To exam- 
ine the diversity and relative abundance of lysine succinylation 
proteins from different growth conditions of Mtb, Western 
blotting analysis was performed using with a polyclonal anti- 
succinyl lysine antibody (Fig. 1). The virulent Mtb strain H37Rv 
were grown in minimal medium, either without an additional 
carbon source or treated with pyruvate, glucose, glycerol, 
succinate, or acetate and harvested at 30 day or 60 day. 
Strong immunoblot signals were observed at both 30 day and 
60 day (Fig. 1), indicating that lysine succinylation was abun- 
dant in H37Rv under different culture conditions. Immunob- 
lots with succinylated or non succinylated BSA confirmed that 
the protein bands detected were specific (Fig. 1). However, 
succinylation signals were stronger when H37Rv were grown 
in medium containing succinate or pyruvate as a sole carbon 
Source than those in no carbon medium (Fig. 1), suggesting 
that different growth mediums may alter the lysine succinyla- 
tion profile. As expected, Mtb exhibited different growth rate 
in media containing different carbon source (supplemental 
Fig. S1). These results implied that the lysine succinylation is 
likely linked with energy metabolism in Mtb. Moreover, as 
shown in supplemental Fig. S2, we obtained a similar result as 
previously reported on lysine acetylation (66) when the strain 
was grown on different carbon source. 

Enrichment and Identification of Lysine Succinylation Sites 
in Mtb—To identify the lysine succinylation sites in Mtb, we 
used a similar procedure as previously described (7, 14, 16) 
(Fig. 2A). In total, we identified 1561 unique modified peptides 
encompassing 1545 succinylation sites (class l) from 626 
proteins with a FDR below 196 for modified peptides. The 
proteomic data showed that 61.496 of identified unique pep- 
tides (2541) across 84.596 of the identified proteins (741) were 
immune-enriched, indicating the highly efficient enrichment 
procedure were used in this study (Fig. 2B and 2C). All spectra 
containing succinylation were manually inspected to ensure 
diagnostic b- or y-ion series as previously described (44, 45). 
Detailed information of all identified succinylated peptides are 
provided in the supplemental Table S1. The raw data has 
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bon sources altered succinylation 
profiles of Mtb H37Rv. Coomassie blue 
staining of protein lysates from A 30 
days and B 60 days of cultures. Coo- 
massie blue staining was used for the 
loading control. Western blotting analy- 
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TABLE | 
Comparison of Mtb H37Rv succinylome with other published succinylomes. “su” = succinylation 


Species Genome size (Mb) Proteins No. of su-sites No. of su-proteins Homologue with H37Rv Ref 
E. coli DH10B 4.47 4,124 2,580 670 128 (14) 
E. coli BW25113 4.64 4,141 2,572 990 153 (16) 
S. cerevisiae S288c 12.16 5,906 1,345 474 60 
H. sapiens (HeLa) 3,209.29 67,778 2,004 738 24 
M. musculus (liver) 2,798.79 73,462 2,140 750 58 
D. melanogaster 139.485 17,526 - 104 6 (7) 
M. musculus (MEFs) 2,798.79 73,462 2,565 779 75 (19) 
M. musculus (liver) 2,798.79 73,462 1,190 252 45 (20) 
M. tuberculosis H37Rv 4.41 4,018 1,545 626 - In this study 


been uploaded onto a public database PeptideAtlas (http:// 
www.peptideatlas.org) and can be accessed through the fol- 
lowing link: http://www.peptideatlas.org/PASS/PASS00518. 
Annotated peptide spectra for all succinylated peptides were 
deposited in PeptideAtlas and supplied as supplementary 
Dataset. The average peptide score was 117.9402 (Fig. 2D) 
and the overall absolute peptide mass accuracy was 0.5752 
parts per million (ppm) (standard deviation, 0.6173 ppm) (Fig. 
2E), confirming the high accuracy of modified peptide data 
obtained from MS. Furthermore, to assess the distribution of 
the average modified peptide length and succinylation sites, 
we calculated the number of amino acids per peptide and 
identified modification sites per protein. Supplemental Fig. 
S3A showed that the average peptide length was 15 amino 
acids and the average site coverage per identified protein was 
three lysine succinylation sites (supplemental Fig. S3B). We 
reasoned that multiple succinylation sites on one protein may 
play an important role in stabilization/destabilization of protein 
conformation and fulfill the delicate regulatory function. 
Conservation of Succinylated Proteins — To explore the ev- 
olutionary conservation of all identified lysine succinylation 
sites in this study, we searched the orthologs of 626 succi- 
nylated proteins of Mtb against the following previously re- 
ported succinylomes: Escherichia coli DH10B (14), Esche- 
richia coli BW25113 (16), Saccharomyces cerevisiae S288c 
(16), Homo sapiens (16), Mus musculus (16), and Drosophila 
melanogaster (7). The orthologs that showed the highest ho- 
mology to each other in both directions were counted and 
reported as the percentage of the total number of proteins in 
all analyzed species. The percentage of Mtb homologs of 
succinylated proteins in diverse organisms is showed in Table 
|. As a result of a global homology search using the amino acid 
sequence of target proteins, orthologs of 205 (—32.796) Mtb 
succinylated proteins can be detected in the succinylomes of 
other species (supplemental Table S2). Moreover, we found 
that 164 and 107 Mtb succinylated proteins have succinylated 
homologs in bacterial and eukaryotic cell, respectively (sup- 
plemental Fig. S4). Interestingly, most succinylated proteins 
are involved in central carbon metabolism and lipids metab- 
olism. These conserved succinylation events suggest that 


succinylation may play an important role in essential biologi- 
cal processes in this pathogen. 

Functional Annotation and Cellular Localization of Succinyl- 
ated Proteins — To better understand the lysine succinylome in 
Mtb, we performed functional annotation analysis using the 
Gene Ontology (GO), Kyoto Encyclopedia of Genes and Ge- 
nomes (KEGG), and Pfam database. Firstly, we investigated 
the GO functional classification of all the succinylated proteins 
based on their biological process and molecular function (Fig. 
3). The classification results for biological process showed 
that the largest protein group of succinylated proteins was 
related to growth (26.6096), followed by cellular process 
(18.4496), metabolic process (17.0496) and response to stim- 
ulus (9.3396). Interestingly, 1.0596 and 0.1296 of the total 
succinylated proteins were involved in immune system pro- 
cess and signaling, respectively, which could be an important 
mechanism for the survival of Mtb. In the molecular function 
classification, most succinylated proteins were related to the 
binding of various targets and enzyme activity. Approximately 
half of the lysine succinylated proteins (50.7696) were cate- 
gorized as catalytic activity-related proteins. Another large 
succinylated protein group is binding to targets (45.45%) 
(supplemental Table S3). Within the cluster of cellular local- 
ization, a total of 74.4496 of the succinylated proteins was in 
the cytoplasmic category by the bacterial protein subcellular 
localization prediction program PSORTb (47). This result was 
also in agreement with the biological process analysis that 
most succinylated proteins were involved in cellular metabo- 
lism. Importantly, we also found that a sizable portion of 
succinylation proteins were located on cytoplasmic mem- 
brane (11.98%), cell wall (0.4896), and extracellular (0.8%). 
These proteins could play a role in regulating the envelope- 
related proteins to impact the virulence of mycobacteria. 

To further elucidate the biological functions of succinylated 
proteins, we performed enrichment analysis using Cytoscape 
plugin BinGO. As shown in supplemental Figs. S5 and S6, the 
succinylated proteins were mostly enriched in growth (p = 
9.86 x 10715), and also statistically enriched in cellular me- 
tabolism process with specific enrichment in cellular ketone 
metabolic process (o = 3.6 x 10 9), carboxylic acid meta- 
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Fic. 3. Histogram representations of the distribution of identified succinylated proteins according to their biological processes, 


molecular functions and cellular localization. 


bolic process (p = 3.86 x 10 9), oxoacid metabolic process 
(o = 3.86 x 10 9), and organic acid metabolic process (p = 
3.86 x 107°) (supplemental Fig. S5 and supplemental Table 
S44). Consistently, the GO enrichment analysis of molecular 
functions further demonstrated that many functions were en- 
riched in our set, including catalytic activity and binding (sup- 
plemental Fig. S6 and supplemental Table S4A). We next 
performed cellular compartment analysis of all succinylated 
proteins and obtained the similar results. As expected, a 
significant portion of succinylated proteins were significantly 
enriched in membrane and cell wall (supplemental Fig. S6 and 
supplemental Table S4A). Moreover, proteins involved in ri- 
bosome, citrate cycle (TCA cycle), fatty acid metabolism, 
glycolysis/gluconeogenesis, pyruvate metabolism, oxidative 
phosphorylation, glyoxylate, and dicarboxylate metabolism 
were also enriched, as determined through the KEGG path- 
way analysis (Supplemental Fig. S7 and supplemental Table 
S4B), indicating that lysine succinylation substrates were as- 
sociated with the ribosome and cellular metabolism events. 
The Pfam domain analysis illustrated that most succinylated 
substrates were enriched with the acyl-CoA dehydrogenase 
domain, which catalyzed the first step in the B-oxidation of 
fatty acids and catabolism of some amino acids (supplemen- 
tal Fig. S7 and supplemental Table S4B). Taken together, 
these findings suggested that the identified succinylated pro- 
teins have a wide distribution of functions and cellular local- 
izations and render particular enrichment to metabolic pro- 
cess in Mtb. 

Sequence Recognition Motifs and Local Structural Proper- 
ties—We further compared the position-specific amino acid 
frequencies of the surrounding sequences (six amino acids to 
both termini) of succinylated lysine residues with those of all 


lysine residues that occurred in the Mtb proteome. We found 
that negatively charged amino acids (glutamic acids) was 
drastically overrepresented at position —1 using motif-X (Fig. 
4A and supplemental Table S5A). Then, we compared our 
data sets to the reported succinylome of E. coli, showing that 
the succinylation pattern found in E. coli BW25113 were quite 
different from Mtb site-specific succinylation motifs. Interest- 
ingly, the same amino acid (glutamic acids) was preferred at 
position +1 in E. coli DH10B. To further assess if there is 
significant enrichment or depletion of specific amino acids 
with respect to the general amino acid composition of the 
entire proteome, a position-specific intensity map was gen- 
erated as previously described (42). The result showed a 
preference for polar acidic and polar basic amino acids (Fig. 
4B). It can be seen that aspartic acid was most commonly 
found at +1 and glutamic acid was most commonly found at 
—1, +2. As well, we observed a significant preference for 
lysine at the —6 to —4, and +5 positions. These results 
indicated that succinylation preferentially occurred in lysine- 
rich regions and polar acidic/basic amino acid regions of 
proteins. We next compared these data with the patterns 
generated from other succinylomes. Although we did not 
detect a strong bias for a specific succinylation site motif in 
E. coli DH10B, there were modest biases for flanking amino 
acid residues. The different preference of amino acid residues 
surrounding lysine sites suggested the unique substrate pref- 
erences in Mtb. Next, we predicted the structural features of 
lysine succinylation sites in identified proteins with NetSurfP 
algorithm (55). As expected, succinylation sites were enriched 
on the protein surface (84.496) as compared with 75.596 of the 
all lysine residue (p = 3.079 x 10 ?). In accordance with 
previous report (16), the succinylated lysine was more fre- 
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Fic. 5. The complete interaction network of identified succinylated proteins and their interacting proteins. The succinylated proteins 
were grouped using functional annotation and interaction network was visualized with Cytoscape. The red node border represents the 


identified succinylated proteins. 


quently found in structured regions. Overall, 54.4% of succi- 
nylated lysines were predicted to be localized in a helix and B 
sheet structures in comparison to 45.5% of all lysines (Fig. 4C 
and supplemental Table S5B). Therefore, lysine succinylation 
may have an effect on the surface properties of modified 
proteins similar as other PTMs, such as acetylation and 
phosphorylation. 

Functional Interaction Networks of Succinylated Proteins in 
Mtb—We generated a protein interaction network of all suc- 
cinylated proteins using protein interaction information. The 
protein interaction map was visualized using Cytoscape (48) 
and consisted of a large network covering 1672 proteins 
among which there were 462 identified succinylated proteins 
(supplemental Table S6). Although further functional experi- 
ments were needed to validate the potential interactions, such 
bioinformatic analysis should be useful for formulating testa- 
ble hypotheses to understand the function of the identified 
succinylated proteins in Mtb. Then, the functional category 
was used to group the identified proteins (supplemental Table 
S7). As shown in Fig. 5, several complexes and cellular func- 
tions formed prominent and highly connected clusters, such 
as cell wall and cell processes, intermediary metabolism and 
respiration, lipid metabolism, and virulence detoxification ad- 


aptation. The interaction networks of succinylated proteins 
involved in the ESX systems as well as virulence, detoxifica- 
tion and adaptation function in Mtb are illustrated in supple- 
mental Fig. S8. We speculate that the physiological interac- 
tion among these protein clusters might contribute to the 
cooperation and/or coordination of their functions in the con- 
trol of diverse signaling and regulatory pathways in Mtb. 

Succinylation of Enzymes Involved in Central Metabolism — 
Increasing evidences show that lysine succinylation may play 
a regulatory role in metabolic pathways, similar with lysine 
acetylation (4, 67). We investigated the succinylation of met- 
abolic enzymes in Mtb by mapping succinylated proteins to 
KEGG metabolic pathways (supplemental Table S8). As 
shown in Fig. 6, a large proportion of metabolic enzymes were 
succinylated. The fact that nearly all enzymes in major meta- 
bolic pathways are subjected to succinylation implies that this 
PTM may regulate cellular metabolic process at multiple lev- 
els in Mtb. For example, a large proportion of metabolic 
enzymes were mapped to two essential pathways for the 
growth and survival of Mtb, i.e. mycolic acid biosynthesis 
(supplemental Fig. S9A) and purine nucleotides de novo bio- 
synthesis pathways (supplemental Fig. S9B), were identified 
to be succinylated. 


Fic. 4. Bioinformational analysis of succinylation sites. A, Sequence Logo representation of significant motifs identified by Motif-X 
software. The motifs with significance of p < 0.000001 are shown. B, Position-specific under- or over-representation of amino acids flanking 
the succinylation sites. Colors were plotted by using intensity map and represent the log,, of the ratio of frequencies within succinyl-13-mers 
versus nonsuccinyl-13-mers (blue shows enrichment, yellow shows depletion). C, Distribution of succinylated and nonsuccinylated lysines in 
protein secondary structures. Probabilities for different secondary structures (o helix, beta-strand and coil) of succinylated lysine were 
compared with the secondary structure probabilities of nonsuccinylated lysine on all proteins identified in this study. Significance was 
calculated by Wilcoxon test. 
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Fic. 6. Central carbon metabolism network in Mtb H37Rv. The identified succinylated proteins were highlighted in red. 


Confirmation of Succinylated Proteins by Immunoprecipita- 
tion and Western blotting — Based on how well their biological 
functions and importance are understood, ten newly identified 
succinylated proteins, including AtpA, AtpF, KasA, DesA1, 
FadD31, Fgd1, GIcB, HspX, MtrA, and Pks2, were selected for 
verification. These proteins were subjected to IP and Western 
analysis (supplemental Fig. S10). The immunoprecipitated 
proteins were first confirmed by their specific antibodies. 
Anti-succinyl lysine antibody detected the respective pro- 
teins, confirming that they contain succinylated lysine resi- 
dues. In comparison, the unmodified BSA (negative control 
protein) was not detected by anti-succinyl lysine antibodies, 
whereas the immuno-reactive signal was strongly observed 
when the succinylated BSA was used as positive control 
(supplemental Fig. S10). These data provided additional bio- 
chemical evidence for these identified succinylated proteins. 

Effect of Succinylation on Acs Activity — Acs plays a crucial 
role in central metabolism and is important for maintaining 
adequate levels of Ac-CoA, a metabolite of many anabolic 


and catabolic processes. It is reported that acetylation of Acs 
completely abolishes the conversion of acetate to Ac-CoA 
(57, 68, 69) (Fig. 7A). In this study, we identified two reliable 
succinylation sites on Acs, Lys193 and Lys366 (Fig. 7B). We 
speculate that the Acs activity may be regulated by reversible 
succinylation. To test our hypothesis, the acs gene of Mtb was 
cloned and was overexpressed as a His-tagged fusion pro- 
tein, which was purified by Ni?* chromatography. As shown in 
Fig. 7C, the recombinant Acs were both acetylated and suc- 
cinylated. Based on the conservation analysis, lysine at 193 
(in coil) and 366 (in o helix) positions of Acs were conserved in 
the Mtb orthologs, suggesting that these residues could be 
important for an evolutionarily conserved function (supple- 
mental Fig. S114). To determine the functional consequences 
of succinylation on Acs acitivity, we converted the modified 
residues at the two positions (Lys193 and Lys366) to arginine 
(R), which essentially locks Acs into a constitutively desucci- 
nylated state. All mutations were verified by both DNA se- 
quencing (supplemental Fig. S11B) and MS analysis 


Molecular & Cellular Proteomics 14.4 


805 


ASBMB 


ow 
pm 
— 


cof 


MOLECULAR & CELLULAR PROTEOMICS 


MCP 


Succinylome Analysis of Mycobacterium tuberculosis 


A 


Glycolysis/gluconeogenesis SDS-PAGE Anti-Acetyllysine 


TCA cycle/glyoxylate bymass 


Mg*/ATP/CoA — AMP/PP, |Lipid synthesis 
Other processes 


Marker Acs BSA ac-BSA Marker Acs BSA ac-BSA 


Acetate Acetyl-CoA 
2*(37)-0- Acetyl-CoA synthetase 


acetyl-ADPr Succinyl-CoA 


Deacetylase 
Desuccinylase 
(CobB) 


Acetyl-transferase 
Succinyl-transferase ? 


Acetyl/Succinyl 
Acetyl-CoA synthetase 


NAD+ 


B SDS-PAGE Anti-Succinyllysine 
Marker Acs BSA su-BSA Marker Acs BSA su-BSA 
localization Es 
G PSPLK 66.993 193 Coil ! 
mvagg;  Acety-CoA K(su) o 
synthetase MFMK(su)WGR 67.035 1 366 Helix 
GK(193, su)PSPLK MFMKW(366, su)GR 


Scan number 2110 Raw file t Method FTMS: HCD Peptide score: 66.99 Sean number 7169 Raw file 1 Method FTMS: HCD Peptide score: 67.03 


st T 


9 


s 
tn > 


Relative Abundance 
0 20 40 60 80100120 
0 sO I SI c sc 
|. 01] Ansuoyug 
Relative Abundance 
0 20 40 60 80 100 120 
t 
I, or] Arsuayug 


100 " 200 " 300 — 400 — 500 — 600 100 ' 300 400 — 500 
m/z [amu] m/z [amu] 


e| 
= 


200 Lanel Lane2 Lane3 Lane4 3.5 4 OLane1 mLane2 
p<0.01 
| en | 


2.5 
p<0.01 


Relative enzyme activity (%) = 
m 
tA 
e 


Relative intensity 


1.5 
T 
0 
Mt-Acs t + + + 0.5 
> y 
FSF SES MtCobB + + + - 
* 5 NADP = = = - i Desuccinylation D lati 
e NAM = M i i esuccinylation Deacetylation 


Fic. 7. Succinylation profiling of Acs of Mtb. A, Schematic illustration of regulation of Acs activity by reversible acetylation and 
succinylation. B, Identification of the modification sites of Acs. The representation of the corresponding MS/MS spectra of succinylated 
peptides from Acs annotated with a comprehensive series of b and y fragment ions. Details of succinylation sites were listed in the 
supplemental Table S1. C, The acetylation/succinylation levels of purified recombinant Acs were determined by immunoblotting using specific 
anti-acetyl/succinyl lysine antibodies. D, Succinylation of Lys193 and Lys366 affect the enzyme activity of Acs. The Acs and its mutants 
(Acs-K193R, Acs-K366R, and Acs-K193R/K366R) were expressed in E. coli BL21 (DE3) and the specific activity of Acs was determined. Data 
are means + S.D. from three independent assays. E, NAD  -dependent deacetylation/desuccinylation of Acs by CobB protein (Rv1151c). Acs 
was incubated with CobB in the presence of NAD* (0.5 mm) for 10h at 25 °C. Nicotinamide (NAM, 2 mm) was added to the reaction as CobB 
inhibitor. F, Quantification of acetylated/succinylated Acs in the reactions from E. Each sample was standardized by comparing to the 
Coomassie blue stained gel. Error bars indicate S.D. of three measurements. 


(supplemental Fig. S11C-S11F). Then, we measured the en- data suggest that succinylation on Lys193 and Lys366 was a 
zymatic activity of Acs and its mutants. The results showed negative regulatory modification on Acs activity. 

that mutation of Lys193 and Lys366 to R (K193R, K366R, NAD* Dependent Desuccinylase Activity of Mtb— The sirtuin 
K193R/K366R), which conserved the positive charge and proteins, which are one conserved family of deacetylases, 
structure with unmodified lysine (15), led to a significant in- have a unique ability to remove acetyl and other acyl groups 
crease of the enzymatic activity (o — 0.01) (Fig. 7D). These from protein lysine side chains in a nicotinamide adenine 
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dinucleotide (NAD*)-dependent manner (70). Recently, the 
human sirtuin protein (Sirt5) was showed to be an efficient 
NAD*-dependent protein lysine desuccinylase and de- 
malonylase (8, 71). In E. coli, a sirtuin-like protein CobB has 
deacetylase/desuccinylase activity (14). In Mtb H37Rv, the 
sirtuin-like protein, CobB (Rv1151c) is able to deacetylate Acs 
(57, 69). Thus, CobB may have desuccinylase activity in Mtb 
H37Rv. To test this possibility, we expressed and purified 
CobB to measure the protein desuccinylase activity in the 
presence of NAD*. As shown in Fig. 7E and 7F, CobB can 
desuccinylate Acs, while no reduction of succinylation was 
observed when CobB or NAD* was withdrawn from the re- 
action. Consistent with previous reports, CobB can also 
deacetylate Acs with similar efficiency (69, 72). Therefore, our 
data suggest that CobB is a bifunctional enzyme with both 
lysine desuccinylation and deacetylation activities in Mtb. 
Molecular Dynamics Simulations of Acs— To gain insight 
into potential mechanisms by which succinylation might affect 
Acs enzymatic activity, we used molecular dynamics (MD) 
simulations to determine how succinylation of K193 and K366 
influences the structure of Acs. It is reported that Acs uses the 
CoA as substrate to form the Ac-CoA complex and the resi- 
dues in Acs that interact with the CoA may affect the binding 
of Acs with CoA by changing the conformation (73). In this 
study, we found the CoA-binding region marked with circular 
section in the unmodified Acs of Mtb, suggesting that the Acs 
could interact with the substrate CoA (Fig. 8A and supple- 
mental Movie S1). When the lysine residues at 193 and 366 
positions were mutated to R, we found that the CoA binding 
region in the mutants and the protein structures were almost 
accordance with the unmodified Acs structure (Fig. 88-8D 
and supplemental Movies S2-S4). Therefore, the mutation of 
lysine residue to R can lock Acs into a constitutively desuc- 
cinylated state and maintain the protein structures to bind 
with CoA. The MD simulations of the succinylated Acs 
(K193su, K366su, K193su/K366su) showed that succinylation 
of K193 and K366 resulted in the change of charge at each 
modified position and led to the conformation change at this 
region, which may prevent the binding of Acs with CoA (Fig. 
8E-8G and supplemental Movies S5-S7). Therefore, succinyl- 
ation at K193 and K366 could change the residue charge 
status from +1 to —1 and consequently disrupt the favorable 
interaction with the near helix region, leading to partially or 
completely perturbation of the neighboring regions and de- 
stabilization of the protein conformation. As shown in Fig. 8H 
and 8/, the structures of succinylated (K193su, K366su, and 
K193su/K366su) and unmodified Acs were significantly differ- 
ent, whereas the R mutants (K193R, K366R, and K193R/ 
K366R) showed nearly identical structures with the unmodi- 
fied Acs during the entire course of simulations. Furthermore, 
K193 and K366, which were positively charged before succi- 
nylation, were located in coil and a helix region, respectively. 
The simulations showed that the mutants and unmodified Acs 
shared the similar secondary structures and underwent little 


change from the starting to their final states. However, the 
secondary structures of succinylated Acs were significantly 
different from that of unmodified Acs (Fig. 8J and 8K). Owing 
to the perturbation of the secondary structures and the con- 
formational fluctuations, the structures of succinylated Acs 
have a strong tendency to form the hydrophobic solvent- 
accessible surface area (Fig. 8L and 8M). Together, these 
observations suggest that the R mutants mimic the structure 
of unmodified Acs at the molecular level and the succinylation 
could influence the conformation of Acs. 


DISCUSSION 


In this study, we report the first extensive data on lysine 
succinylome in the virulent Mtb strain H37Rv. A total of 1561 
unique modified peptides encompassing 1545 succinylation 
sites (class l) were identified from 626 H37Rv proteins. Be- 
cause lysine succinylation has emerged as an important reg- 
ulatory PTM and may impact diverse metabolic pathways, our 
data set provides a rich resource of putative regulatory mod- 
ification sites in this life-threatening pathogen. We show that 
proteins in major metabolic pathways are targeted by succi- 
nylation. These pathways include glycolysis, fatty acid oxida- 
tion, TCA cycle, and mycolic acid biosynthesis, which are 
crucial components of Mtb's energy metabolic networks. 

Among all metabolic pathways, mycolic acid biosynthesis 
pathway is unique to mycobacteria and is crucial for their 
success as pathogens (74). Our data showed that a large 
proportion of enzymes in mycolic acids biosynthesis pathway 
were succinylated (supplemental Fig. S9A). For example, the 
single fatty acid synthase (FAS) | type protein FAS, which 
generates CoA esters from Ac-CoA primers, creating precur- 
sors for elongation by all other fatty acid and polyketide 
systems, was found to be succinylated. In FAS-II system of 
mycolic acid biosynthesis, sets of enzymes, which can be 
inhibited by many drugs used in treatment of TB (75), were 
also succinylated. Therefore, the reversible succinylation, al- 
though the mechanism remains unknown, could be an impor- 
tant mechanism for regulating the interaction of these drugs 
with FAS-ll-related enzymes, and could also be targets for 
novel drug therapies. The recent observation showed that 
Mtb is metabolically flexible and can even alter the host 
cellular metabolism, such as the pathogen-induced dysregu- 
lation of host lipid synthesis, to create a niche that is ideally 
suited to its persistent lifestyle (76). Mtb can utilize the lipid 
bodies accumulated in the infected host as a nutrient source 
via B-oxidation pathway and this seems to sustain intracellular 
Mtb in a physiological state. Thus, aside from fatty acid bio- 
synthesis pathway, the f-oxidation of fatty acids, which 
Serves as carbon and energy source, is also required for 
survival of Mtb in the infected host cell (77). Several enzymes 
involved in this pathway were found to be succinylated in this 
study (Fig. 6), indicating that lysine succinylation may play a 
crucial role in regulating the fatty acid metabolism of Mtb 
during infection. 
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Fic. 8. Structural modeling of Acs. Final structures of the simulations of Acs were shown. The circular section marks the CoA-binding 
region A-D. A, Nonsuccinylated (Wt) Acs. Site-directed mutants of B, Acs-K193R; C, Acs-K366R; and D, Acs-K193R/K366R. Acs succinylated 
at E, K193; F, K366; and G, K193 and K366. H-/, Analysis of the root mean square deviations (RMSDs), root mean square fluctuations (RMSFs) 
and radius of gyrations (Rgs) for the different Acs systems. (Top) Trajectories of the overall RMSDs of the different Acs systems. (Middle) Radius 
of gyrations profiles of the different Acs systems. (Bottom) Residue-specific RMSF profiles of the different Acs systems. J-K, Secondary 
structure evolution of Acs and its mutants. The evolution in secondary structure at each frame was monitored using the dictionary of protein 
secondary structure (DSSP) algorithm. In the stripes each pixel represents the secondary structure (color-coded) of a residue (185-190 or 
355-370, x-dimension) at a given time in simulation (y-dimension). L-M, Time evolution of the solvent accessible surface area (SAS) calculated 


for the different systems that were simulated. 


It is now widely accepted that metabolic adaptation to the 
host environment is a defining feature of the pathogenicity of 
Mtb (38, 78). However, we lack biochemical knowledge of its 
metabolic networks and how they are regulated. The succinyl- 
CoA is a succinyl donor to lysine and the concentration of 
succinyl-CoA may regulate succinylation levels globally. Our 
understanding of the role of succinyl-CoA and lysine succinyl- 
ation on the control of Mtb metabolism remains rudimentary. 
However, their importance in the control of metabolism in bac- 
teria, yeast and human, and the conservation of the enzymes 


involved in these pathways in Mtb, suggest that lysine succinyl- 
ation may work together with many other PTMs to regulate the 
metabolic pathways of Mtb. Mtb may have unique mechanisms 
involved in controlling metabolic flux via succinyl-CoA and pro- 
tein succinylation, because it resides in humans as its only 
known natural host and a large part of the energy status in Mtb 
occurs through the interaction with the human cells (35). Re- 
versible succinylation is a strong candidate for dynamically con- 
trolling metabolic pathways and enabling Mtb to adapt to the 
changing host environment. To better understand the role of 
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succinyl-CoA and protein succinylation in Mtb metabolism, it is 
essential to acquire an in-depth description of the succinylome 
in Mtb. Furthermore, it is important to understand how the 
succinylation of metabolic enzymes regulates metabolic path- 
ways in Mtb. More studies need to be done to accomplish this. 
Firstly, examine how environmental conditions as well as the 
carbon and energy balance of Mtb cells can influence the level 
of lysine succinylation. Secondly, determine the effects of suc- 
cinylation on enzymatic activity. Thirdly, identify lysine succinyl- 
transferase and desuccinylase involved in controlling the suc- 
cinylation status of proteins in Mtb. In mammalian cell, the 
histone H3 propionylation might be generated by the histone 
acetyltransferase p300 and that selection of donor molecules 
(propionyl-CoA versus acetyl-CoA) may determine the differ- 
ence of modifications (79). In this study, we identified the first 
desuccinylase, CobB, in Mtb, however, it remains unknown 
whether succinyltransferase and more desuccinylases exist in 
Mtb. Acetylation of Acs by the acetyltransferase (Pat) inacti- 
vated the enzyme in S. enterica, whereas deacetylation by the 
CobB reactivated it (80). We found the Pat homolog in Mtb 
(Rv0998) and further studies are thus necessary to identify the 
succinyltransferase in Mtb. 

Being only one of many PTMs, lysine succinylation com- 
petes with other modifications, such as acetylation, for the 
same lysine residues. Accumulating evidence indicated that 
lysine acetylation is a conserved PTM that links Ac-CoA me- 
tabolism and cellular signaling (81, 82). Acs is the enzyme 
catalyzes the conversion of acetate into Ac-CoA and the level 
of Ac-CoA can fluctuate according to the changing metabolic 
conditions (83). In this study, we demonstrated that the en- 
zymatic activity of Acs is affected by the lysine succinylation. 
It is thus expected that crosstalk exists between protein suc- 
cinylation and acetylation may modulate protein activity and 
metabolic pathways in both Mtb and human cells. Jungblut 
and Schlüter proposed a systematic nomenclature for the 
comprehensive description of protein species and suggested 
that protein species sequence including PTMs should replace 
the basic unit protein sequence (84, 85). Therefore, our suc- 
cinylome data may help to provide the accurate and detailed 
description of protein species in Mtb and to elucidate poten- 
tial biological roles of each protein species. 

In summary, we report the first large-scale, high resolution 
MS-based survey of lysine succinylation in Mtb. The provided 
dataset may serve as important starting point for assessing 
the functional relevance of succinylation sites in this life- 
threatening pathogen. Given that the high abundance of lysine 
succinylation and its induced conformation fluctuation, it is 
likely that reversible succinylation may be an important regu- 
latory mechanism in metabolism. It would be interesting and 
important to investigate the biological functions of reversible 
succinylation in metabolic regulation in Mtb. 
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